Measurement of hadronic cross section and preliminary results on the 
pion form factor using the radiative return at DA$NE 

The KLOE Collaboration * 
presented by G. Venanzoni a 

a Universita and Sezione INFN Pisa 

In the fixed energy environment of the e + e~ collider DA<I>NE, KLOE can measure the cross section of the 
process e + e~ — > hadrons as a function of the hadronic system energy using the radiative return. At energies 
below 1 GeV, e + e~ — > p — > -k + tt~ is the dominating hadronic process. We report here on the status of the 
analysis for the e + e~ — > n + n~"/ channel, which allows to obtain a preliminary measurement of the pion form 
factor using an integrated luminosity of ~ 73p6 _1 . 



1. Measuring the hadronic cross section 
with the KLOE detector 

1.1. Motivation 

The measurement of the hadronic cross section 
at low energy is of great importance for the im- 
provement of the theoretical error of the anoma- 
lous magnetic moment of the muon, = (g^ — 
2)/2. The hadronic contribution a^ adr is given 
by the hadronic vacuum polarization and cannot 
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be calculated at low energy in the framework of 
perturbative QCD. Following a phenomcnological 
approach, the hadronic contribution can however 
be evaluated from the measurement of R through 
a dispersion relation. 

a H« d r = ( ^)2 r rfs ^)f( g ) ; (1) 
-i 7-)/ \ fT(e"*"e — — rhadrons) jj.ii l 

where R(s) = — 4Taij(3) 1 and the kernel 

3s 

K(s) is a smooth bounded function growing from 
0.63 at threshold to 1 at oo. Due to the 1/s 2 
dependence in the integral, hadronic data at 
low energies are strongly enhanced in the con- 
tribution to a!}°- dr . The error of the hadronic 
contribution is therefore given by the limited 
knowledge of hadronic cross section data. This 
error is the dominating contribution to the total 
error of a^ eo (Saf eo « 5a* adr ) [Q 

The discrepancy between the theoretical calcu- 
lation , and the recent experimental value 
of [p[ depends if r data are used or not in 
the evaluation of a h ^ d ; a recent re-evaluation 
of the hadronic contribution to Q found: 
a theo_ a ex P = (33.9 ±n.2) x 1(T 10 [e+e- based], 
a theo_ a ex P = (16.7 ± io.7) x 1CT 10 [t based], cor- 
responding to 3.0 and 1.6 standard deviations re- 
spectively, reflecting also a 1.6 a discrepancy be- 
tween e + e~ and r evaluations of a^ ad . A discus- 
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sion on this subject can also be found in |9j]. 

This unclear situation makes a precise mea- 
surement of the hadronic cross section at low 
energy mandatory, particularly in the channel 
e + e~ — ► p — » 7r + 7r~ , which is the main ingredient 
for a^ ad . Such a cross section has been already 
measured by the CMD-2 collaboration with an 
energy scan in the range below 1.4 GeV Jh}; an 
accuracy of 0.6% was recently achieved in the re- 
gion between 0.61 and 0.96 GeV 0. 

1.2. Radiative Return 

DA$NE is an e + e~ storage-ring collider 
working at the </> resonance (1020 MeV). As an 
experiment at a collider with a fixed centre of 
mass energy, KLOE can measure the hadronic 
cross section a(e + e~ — > hadrons) as a function 
of the hadronic system energy using the radiative 
return Jl3j , |l4| ] i. e. studying the process e + e~ — > 
hadrons+7. The emission of one photon before 
the beams interact (Initial State Radiation, ISR 
in the following) lowers the interaction energy and 
makes possible to produce the hadronic system 
with an invariant mass varying from the <j) mass 
down to the production threshold. 

The method represents an alternative approach 
to the conventional energy scan used so far for 
hadronic cross section measurements. A very 
solid theoretical understanding of Initial State 
Radiation (described by the radiation function 
H) is mandatory in order to extract the cross 
section er(e + e~ — > hadrons) as a function of 
Qhad f rom the measured differential cross section 
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where Qf lad and # 7 are respectively the invariant 
mass squared of the hadronic system and the ac- 
ceptance cut on the polar angle of the photon. 

Radiative corrections for e + e~ — > 7r + 7r~7, have 
been calculated up to ISR NLO by different the- 
oretical groups IHJllQIIJlll^], and recently 
implemented in a new Monte Carlo generator 
named PHOKHARA [plpl, This generator with 



a claimed accuracy of 0.5% has been used in the 
present analysis. 

We want to stress that the radiative return has 
the merit compared with the conventional energy 
scan that the systematics of the measurement 
(e.g. normalization, beam energy) are the same 
for any experimental point and must not be evalu- 
ated at each energy step. However a precise deter- 
mination of the angle of the hard photon as well as 
the full control of events with the photon emitted 
in the final state (pure QED FSR or other reso- 
nant processes, as for example <p — > 7r + 7r~7 p3[), 
is required. 



2. Analysis of the process 



7T+7T 7 



Due to the importance of the 27r final state for 
a£ ad , we concentrated at KLOE on the analysis 
of the radiative process e + e~ — > 7r + 7r~7, whose 
feasibility has been extensively studied by Monte 
Carlo Q. 

2.1. The KLOE Detector 

KLOE |2j| is a typical e + e~ multiple pur- 
pose detector with cylindrical geometry, consist- 
ing of a large helium based drift chamber (DC, 
p6fl ) , surrounded by an electromagnetic calorime- 
ter (EmC, and a superconducting magnet 
(B = 0.52 T). The detector has been designed 
for the measurement of CP violation in the neu- 
tral kaon system, i.e. for precise detection of the 
decay products of K$ and Kl- These are low 
momenta charged tracks (ir^ , //^ , with a mo- 
mentum range from 150 MeV/c to 270 MeV/c) 
and low energy photons (down to 20 MeV). 
The DC dimensions (3.3 m length, 2 m radius), 
the drift cell shapes (2x2 cm 2 cells for the inner 12 
layers, 3x3 cm 2 cells for the outer 46 layers) and 
the choice of the gas mixture (90% Helium, 10% 
Isobutane; X — 900 m) had to be optimized for 
the requirements prevailing at a factory. The 
KLOE design results in a very good momentum 
resolution: cr p± /p± < 0.4% at high tracking effi- 
ciencies (> 99%). 

The EmC is made of a matrix of scintil- 
lating fibres embedded in lead, which guar- 
antees a good energy resolution <je/E = 
5.7%/ \J E(GeV) and excellent timing resolution 
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Figure 1. KLOE event display of a 7r + 7r 7 event. 



a t = 57ps/^E(GcV) © 50 ps. The EmC consists 
of a barrel and two endcaps which are surround- 
ing the cylindrical DC; this gives a hermetic cov- 
erage of the solid angle (98%). However, the ac- 
ceptance of the EmC below « 20° is reduced due 
to the presence of quadrupole magnets close to 
the interaction point and does not allow to mea- 
sure e.g. the photon of 7r + 7r~7 events with low 
9 j angles. 

It will be shown in the following that an ef- 
ficient selection of the ir + ir~j signal is possi- 
ble, without requiring an explicit photon detec- 
tion. The relatively simple signature of the signal 
(2 high momentum tracks from the interaction 
point, see Fig. [T]) and the good momentum reso- 
lution of the KLOE tracking detector allow us to 
perform such a selection. 

2.2. Signal selection 

The selection of e + e~ — > tt + tt~j events is done 
in the following steps: 

• detection of two charged tracks, with polar 
angle bewteen 40° and 140°, coming from 
a vertex in the fiducial volume R < 8 cm, 
\z\ < 15 cm. The cuts on the transverse 
momentum px > 200 MeV or on the lon- 



gitudinal momentum \p z \ > 90 MeV reject 
tracks spiralizing along the beam line, en- 
suring good reconstruction conditions. The 
probability to reconstruct a vertex in the 
drift chamber is ~ 95% and has been stud- 
ied with Bhabha data, selected using the 
calorimeter only. The overall tracking re- 
construction efficiency has been also evalu- 
ated using n + ir~iT Q events selected by de- 
tecting 7T° in the electromagnetic calorime- 
ter; 

• identification of pion tracks: a Likelihood 
Method, using the time of flight of the par- 
ticle and the shape of the energy deposit in 
the electromagnetic calorimeter, has been 
developed to reject the e + e~ — > e + e~7 
background. A control sample of 7r + 7r~7r° 
has been used to study the behaviour of 
pions in the electromagnetic calorimeter, 
since their interaction are not well repro- 
duced by Monte Carlo. The effect of the se- 
lection on 7r + 7r~7 and e + e~7 events is vis- 
ible in the upper plot of Fig. ||, as a func- 
tion of the track mass, Mtrack'- this vari- 
able is calculated from the reconstructed 
momenta, p+, p-, applying 4- momentum 
conservation, under the hypothesis that the 
final state consists of two particles with 
the same mass and one photon. For each 
event class (/.t + /i~7, 7r + 7r~7) the track mass 
distribution is peaked at the proper mass; 
7r + 7T~7r events populate the region on the 
right of the 7r + 7r~7 peak. As it can be noted 
from Fig. ^, up, e + e~7 events are drasti- 
cally reduced, while the 7r + 7r~7 peak is es- 
sentially unaffected by this selection (the 
estimated signal efficiency is larger than 
98%). 

• cut on the track mass: events are 
rejected by a cut at 120 MeV in the track 
mass. The discrimination between tt and \l 
using calorimeter information is not help- 
ful since pions behave frequently like min- 
imum ionizing particles. After this cut we 
find a contamination of /^ + /i~7 background 
smaller than 1%. 7r + 7r~7r events are re- 
jected with a cut in the two-dimensional dis- 
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Figure 2. Up: track mass distribution before and 
after the likelihood selection. /i + y«~7, Tr + ir~j dis- 
tributions are peaked at the proper mass, three 
pion events populate the region on the right of 
the 7r + 7r~7 peak. Down: cut in the plane (Q^, 
track mass) used to reject 7r + 7r~7r° events. 



tribution of the track mass versus the two 
pion invariant mass squared, Q\^, shown 
in Fig. |[ down. Due to the large produc- 
tion cross section for these events (BR(<fi — > 
tt+tt-tt ) - 15%, a ~ 500 nb) and the 
similar kinematics, some residual contam- 
ination is expected at small values 
{Ql n < 0.4 GeF 2 ). The efficiency of the 
track mass cut, as evaluated from Monte 
Carlo, is - 90%; 

• definition of the angular acceptance: The 
polar angle of the photon, <9 7 is calculated 
from the charged tracks as 180°— 9^, where 
9 n n is the polar angle of the two pion sys- 
tem. Two fiducial volumes for the analy- 
sis of 7r + 7r~7 events are defined: 9 1 < 15° 
or 9 1 > 165° (small angle analysis) and 
55° < 9 1 < 125° (large angle analysis). 
These two regions mainly differ for FSR and 
background contamination p8| ]. 

In this paper we will show the preliminary results 
obtained on the analysis at small angle (# 7 < 15° 
or 9 1 > 165°, 40° < 6* < 140°). 

3. Final State Radiation 

It is now a common understanding that the ra- 
diation of one hard photon by final pions should 
be included in the measurement of the hadronic 



cross section for the evaluation of 



9U33C 



-,had 

Duutiuii lkjl tile cvaiuatiun Ul uL ^ 

However JL4J, in the case of radiative re- 
turn, the differential cross section for the process 
e + e~ — > 7r + 7r~7 is proportional to the pion form 
factor squared at the measured Q 2 ^ only for ISR: 

dcr 



\FAQ 



2\\2 



(<n (2) 



dQ 2 / ISR 

while for FSR the differential cross section Q is 
proportional to the pion form factor evaluated at 
the <j) mass squared: 



/ dcr 



\FM)\^a e 



(Q 2 )-(3) 



\dQ 2 j FSR 

For this reason events due to FSR must be 
properly rejected by the analysis cuts. 

2 If no additional photons are emitted in the initial state. 
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KLOE 7C7ty events after selection -73pb~ 




Figure 3. Ratio FSR/{FSR + ISR) as a function 
of the two pions invariant mass for 6> 7 < 15° or 
1 > 165°, 40° < 6* < 140°. 
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Fig. | shows the ratio FSR/(ISR + FSR) as 
a function of the Q\ v for the small angle region 
(0 7 < 15° or 9 1 > 165°, 40° < 6„ < 140°). As 
it can be noted the average value of this ratio is 
below 0.5%. 

4. Results 

4.1. Effective cross section 

Fig. |J up shows the number of irirj events se- 
lected inside the small angle acceptance cuts at 
the end of the selection chain, for 100 bins be- 
tween 0.02 and 1.02 GeV 2 . This distribution cor- 
responds to ~ 73p& -1 of analyzed data out of 
~ 200 pb^ 1 collected in 2001 and already recon- 
structed. More than 1083000 were selected, cor- 
responding to 15000 events/p& -1 . 

Even before unfolding the spectrum for the de- 
tector resolution effects, the p — oj interference 
(as well as the radiative tail) can be clearly seen, 
showing the excellent momentum resolution of 
the KLOE DC. 

4.2. Pion form factor calculation 

Neglecting FSR interference, the pion form fac- 
tor can be extracted as a function of Q^, from 
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Figure 4. Up: Number of ir + ir 7 events selected 
by the analysis in the region 7 < 15° or B 1 > 
165°, 40° < 0„ < 140°. Down: Global efficiency 
as a function of Q\^- 
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the observed it + it 7 spectrum (N obs ) as: 
N obs 

\ F ^\ 2 = e(Q 2 )LH(Q 2 ) (4) 

where e(Q 2 ) is the global analysis efficiency, L is 
the integrated luminosity and H(Q 2 ) is the NLO 
cross section for e + e~ — > ir + ir~j (only ISR) un- 
der the assumption of pointlike pions. 

These quantities were evaluated in the follow- 
ing ways: 

1. global efficiency: this is the product of 
trigger, reconstruction, filtering, likelihood 
and Mtrack cut efficiency. Apart from the 
latter which was evaluated only by Monte 
Carlo, the other efficiencies were evaluated 
from unbiased samples of data with simi- 
lar kinematics (like 7r + 7r~7r°, Bhabha's or 
7r + 7r~7 itself). The Monte Carlo was used 
to generate the appropriate 7r + 7r~7 phase 
space. 

Fig. [|, down, shows the dependence of 
the global efficiency as a function of re- 
constructed Q 2 ^ (only statistical errors are 
shown). Above OAGeV 2 the average er- 
ror is 2%, mainly dominated by the limited 
Monte Carlo statistics. Below OAGeV 2 the 
error is larger, reflecting again the limited 
Monte Carlo statistics used in the evalua- 
tion of M^ack cut efficiency 

2. Luminosity: The DA$NE accelerator 
does not have luminosity monitors at small 
angle due to the existence of focusing 
quadrupole magnets very close to the in- 
teraction point. The luminosity is there- 
fore measured using Large Angle Bhabhas 
(55° < 6+ _ < 125°, a = 425nb). The num- 
ber of LAB candidates are counted and nor- 
malized to the effective Bhabha cross sec- 
tion obtained from Monte Carlo. 

The precision of this measurement depends 
both on the understanding of experimental 
efficiencies and acceptances and on the the- 
oretical knowledge of the process. 

The systematic errors arising from the LAB 
selection cuts are well below 1%. All the se- 



lection efficiencies concerning the LAB mea- 
surement (Trigger, EmC clusters, DC track- 
ing) are above 98% and well reproduced by 
the detector simulation. The background 
due to /i + //~7, 7r + 7r~7 and 7r + 7r~7r is be- 
low 1%. 

KLOE uses two independent Bhabha event 
generators (the Berends/Kleiss genera- 
tor |H modified for DA<f>NE || and 
BABAYAGA ||). 

The very good agreement of the experi- 
mental distributions (6+.-, with the 
event generators and a cross check with an 
independent luminosity counter based on 
e + e~ — > 77(7) indicates a precision of bet- 
ter than 1%. 

More systematics checks (e.g. the effect of 
a varying beam energy and of a displaced 
beam interaction point) are under way. 

3. determination of H: The NLO cross sec- 
tion for pointlike pions H was computed for 
each bin of Q 2 ^ using a modified version of 
PHOKHARA, in which F„(Q 2 ) = 1. 

2 • 10 6 events were generated within the ac- 
ceptance of the small angle analysis: 6> 7 < 
15° or 6 1 > 165°, 40° < 6 W < 140°. Fig. | 
shows a comparison between the differen- 
tial cross section generated by Monte Carlo 
with a given parametrization of the pion 
form factor (white dots), and the same as- 
suming pointlike pions (black dots). 

The systematical error on the numerical 
evaluation of H for the given acceptance 
cuts is estimated to be less than 1% by 
comparing the pion form factor obtained 
by dividing bin-by-bin the two distribu- 
tions of Fig. H with the expected analytical 
parametrization. 

The pion form factor extracted using eq. |], is 
shown in Fig. |(| The spectrum can be divided in 
three regions according to the measurement error: 
below OAGeV 2 with an error of 5-10%, between 
OAGeV 2 and 0.5 GeV 2 , with an average error of 
3%, and above 0.5 GeV 2 with an average error of 
2%. The measurement error is dominated by the 
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Figure 5. Differential tt + tt~j NLO cross sec- 
tion from Monte Carlo, for a given F 7T (Q 2 ) 
parametrization {white dots), and assuming 
pointlike pions, i.e. F 7r (Q 2 ) — 1 (black dots). The 
cuts are 7 < 15° or 6 1 > 165°, 40° < 9 n < 140°. 



limited Monte Carlo statistics used in the evalu- 
ation of the efficiencies. 

4.3. Fit of the pion form factor 

A preliminary fit was applied to data using the 
parametrization for F V (Q 2 ) found in |33] (KS): 



FAQ 2 



BW t 



l + Q 



(5BW. 



p' 



1 + 



(5) 



The mass and the width of the p as well as a and 
were free parameters of the fit, while the other 
paramaters were kept fixed (to the values of jnj ) . 

The results obtained from the fit were: M p — 
(772.6 ± 0.5) MeV, T p = (143.7 ± 0.7) MeV, 
a = (1.48 ± 0.12) ■ 10- 3 , (3 = -0.147 ± 0.002. 
These values, even if preliminary, are in good 
agreement with the ones found by CMD-2 |[l| 
using a light different parametrization (Gounaris- 
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Figure 6. Up: Preliminary measurement of the 
pion form factor. Overimposed are the analytical 
\F W (Q 2 )\ 2 parametrizations, obtained by a fit to 
the KLOE data using KS model (blue line), and 
the parametrization found from CMD-2 with the 
GS model (red line). The two curves are very 
similar. Down: zoom of \F 7r (Q 2 )\ 2 around the p 
peak. The blue line is obtained by a fit to the 
KLOE data using KS model, while the red one is 
the parametrization found from CMD-2. 
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Sakurai). This agreement can be also confirmed 
from Fig. ||, which shows the results of the two 
parametrization (KLOE blue line, CMD-2 red 
line) overimposed to our experimental data. The 
discrepancy at low could be due either to 
the model inadequacy or to residual background 
events. 

5. Conclusion and outlook 

The status of e + e~ — > 7r + 7r~7 analysis and the 
preliminary results obtained on the pion form fac- 
tor were presented. The analysis was based on an 
integrated luminosity of ~ 73 pb^ 1 , which corre- 
sponds to about 1 /3 of the data collected (and al- 
ready reconstructed) in 2001. The average error 
for single point above OAGeV 2 is ~ 2%, mainly 
dominated by the the limited Monte Carlo statis- 
tics in the efficiency evaluation. FSR and other 
background channels (as for example Bhabha's 
or 7r + 7r~7r°) are kept below 1% by a likelihood 
method and kinematical cuts of the analysis. 
Other sources of systematical errors due to lu- 
minosity, radiative corrections and resolution are 
still under study, and should be included into our 
errors. 

A comparison of the pion form factor with the 
fitted parametrization obtained by the CMD-2 
collaboration, shows a good (preliminary) agree- 
ment, even without unfolding the spectrum for 
the detector resolution effects. Such a compar- 
ison, even at an accuracy of 2-3%, has become 
more and more important in view of the 1.6 a dis- 
crepancy between e + e~ and r data in the recent 
evaluation of a\ ad [|. 

In order to improve the accuracy on a ! ^ ad a fi- 
nal precision for this measurement below 1% is 
needed. This is very important both in the en- 
ergy region around the p, where the CMD-2 col- 
laboration has recently reached a 0.6% accuracy 
as well as below 0.6 GeV, which contributes to 
with about 100 x 10 -10 and is known with a 
worse accuracy |hJ. KLOE can study both these 
regions using 7r + 7r~7 events with the hard photon 
emitted at small and large angle. However in or- 
der to reach such an highly demanding accuracy 
many factors must be well under control, both ex- 
perimentally and theoretically. The results here 



presented, together with the larger statistics up 
to now collected by KLOE (~ 500 pb- 1 ) and the 
intense theoretical work on radiative corrections 
from different groups, are a promising indication 
for achieving such a challenging task. 
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